THE GEODESICS OF ROLLING BALL SYSTEMS 
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Abstract. A well-known and interesting family of sub-Riemannian space are the systems 
involving two balls rolling against each other without slipping or twisting. In this note, we 
show how the sub-Riemannian geodesies of these space, when the two balls are embedded 
bJQf in R 3 xl 3 , are horizontal curves on the intersections of these balls with Euclidean 5-planes. 
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1. Introduction 

Rolling ball systems, in which two spheres with different radii roll against each other 
without slipping or twisting, are particularly interesting examples of sub-Riemannian man- 
ifolds. E. Cartan in |CarlO] was the first to study these configuration spaces, albeit in a 
different form. After the re-emergence of sub-Riemannian geometry as a topic of interest in 
the late twentieth century, rolling ball systems and sub-Riemannian manifolds with growth 
vector (2,3,5) were studied by R. Bryant and L. Hsu in [BH93j. R. Montgomery's invalu- 
able monograph on sub-Riemannian geodesies, |Mon02] . included a section on rolling ball 
systems, and an open question posed by him in that section about the relationship between 
the symmetries of rolling ball systems and the exceptional Lie group G 2 sparked a great deal 
of research over the past decade. See |Zel06| . |Agr07| , |BM09] . and most recently |BH12] for 
more details. 

Surprisingly, with all of this interest in rolling ball systems, we have yet to see an explicit 
description of the sub-Riemannian geodesies of a rolling ball system beyond the singular 
geodesies first discovered by A. Agrachev and A. Sarychev in [AS95j . It turns out that, 
like the Heisenberg group, the sub-Riemannian geodesies of a rolling ball system have a 
very simple characterization. Specifically, if we consider each sphere as embedded in M 3 and 
centered at the origin, and the point of contact between the two spheres as an ordered pair 
111 I 3 x K 3 , then for any given sub-Riemmanian geodesic of our rolling ball system, there 
exists a 5-plane in M 3 x M 3 such that all points of contact on the geodesic lie on this 5-plane. 
Since the rolling ball system comes equipped with a two-dimensional horizontal distribution, 
for almost all 5-planes passing through a point, a unique horizontal curve is determined by 
the 5-plane. Thus this characterization of the sub-Riemannian geodesies of our system is 
quite specific. 

The proof below is quite simple, and should be understandable for an advanced under- 
graduate student. Essentially, we choose five arbitrary points of contact along our candidate 
geodesic 7 and bend the curve at these points. We then show that the coordinates of those 



Date: August 27, 2012. 

2000 Mathematics Subject Classification. Primary: 53C17; Secondary: 53C22. 

1 



2 



DANIEL R. COLE 



points of contact in R 3 x M 3 must be linearly dependent in order for 7 to be a length min- 
imizing geodesic. We note that a similar proof, using only three points, can be used to 
characterize the sub-Riemmanian geodesies of the Heisenberg group (and thus offers yet 
another proof of the isoperimetric problem in the plane). 

2. Results 

To begin, let Mi and M 2 be the two spheres in our rolling ball system. Without loss of 
generality, we may assume that M 1 has radius 1. Let r be the radius of M 2 . Each sphere sits 
in its own copy of Euclidean space centered at the origin. Therefore every point in Mi can 
be represented by a unit vector u, and every point in M 2 can be represented by rv, a unit 
vector multiplied by r. A point of contact between Mi and M 2 consists of an ordered pair 
(u, rv) of points in Mi and M 2 , respectively. In general, a state in our rolling ball system 
can be represented by a point of contact and two unit vectors in M 3 , a and b, where u • a = 
and v • b = 0. We can think of a as a unit tangent vector to Mi at u, and b as a unit tangent 
vector to M 2 at rv. If we map Mi and M 2 into a common copy of M 3 so that the images of 
Mi and M 2 are tangent to each other at the images of u and rv, then the images of a and 
b coincide with each other on the common tangent plane. Of course, we can choose another 
pair of unit tangent vectors to represent the same state: u x a and — v x b work nicely 
(if we assume that the two sphere images lie in each other's exteriors, as opposed to one 
sphere image nested inside the other). For each point of contact (u, rv), we can define an 
equivalence relation: two pairs of unit tangent vectors (ai,bi) and (a 2 ,b 2 ) define the same 
state at (u, rv) if and only if both (a l5 a 2 ) = (b 1; b 2 ) and (ai, u x a 2 ) = — (bi, v x b 2 ). We 
can define M to be the six dimensional manifold consisting on all 4-tuples (u, v, a, b), and 
we can define M, our rolling ball system, to be the five dimensional manifold consisting of 
all equivalence classes, or states, of those 4-tuples. 

The manifold M inherits a natural sub-Riemmanian structure with growth vector (2, 3, 5) 
from the Riemannian geometry of the spheres. We interpret the "no slipping, no twisting" 
restriction on the motion of the two sphere to mean that a horizontal path 7 in M is 
completely determined by motion along either sphere. Specifically, for a 4-tuple (u, v,a, b), 
a unit speed path along the great circle on Mi in the direction of a cos 9+ (u x a) sin 9 will force 
unit speed motion along the great circle on M 2 in the direction of b cos 9 — (v x b) sin 9. This 
motion will include the parallel transport of a and b induced by the Levi-Civita connections 
on Mi and M 2 . Altogether, a horizontal direction on M will be a linear combination of 

the vectors Xi = ^a, b, — u, j and X 2 = (u x a, — v x b, 0, 0), which together form an 

orthonormal frame. The distribution formed by these two vector fields is compatible with 
the equivalence relation defined above and therefore defines a two-dimensional horizontal 
distribution on M. 

One of the nice features of M is that it inherits a great deal of symmetry from Mi 
and M 2 . Specifically, if gi,g 2 G 5*03, then the map (u, v,a, b) 1— > (giU, g 2 v, gi&, g 2 h) is 
an isomorphism on M, and therefore on M. With this observation in mind, it suffices to 
describe the horizontal geodesies from any state p = [(u , v , a , b )] in M to the state 



n= [((0,0,1), (0,0,1), (1,0,0), (1,0,0))] 
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the state where the spheres are in contact at their respective north poles, and the unit 
tangent vectors at these north poles both point towards the positive x-axis. We can define 
a coordinate system (wi, W2, w 3 , W4, W5) for M centered at n as follows: if (u, v,a, b) is a 
representative of a state near q, then w\ — ui, W2 — U2, w 3 — v±, — V2, and 



It is easy to show that the above coordinates do not depend on the representative used for 
the state. 

Let j(t) be an absolutely continuous horizontal path on M with a unit speed parametriza- 
tion such that 7(0) = p and there exists T such that 7(T) = n. Assume that 7(f) extends 
past n to some point q, and that j(t) is a length minimizing geodesic from p to q. 

Arbitrarily fix five times ii, £2, £3, £4, £5 on (0,T). We are going to perturb 7 by bending 
its image at each of these five times. Specifically, let R 6 U G S0 3 be the counterclockwise 
rotation about the axis passing through the unit vector u by an angle 9 (here, the rotation 
is counterclockwise as viewed with u pointing towards us). Let a = (cui, a 2 , ct 3 , «4, 0:5) be a 
5-tuple of angles. Define the the path -y a (t) = [(u Q (t), v Q (f), a Q (t), b a (t))] as follows 



w 5 = (01U3 - a 3 ui) b 2 + (61 u 3 - b 3 U2) a 2 



< 



u(t) 

i? u(t 1 ) i? u(t 2 ) i? u(t 3 ) U ( t ) 

R% ,R% ,R a % ,R a % .u(t) 

u(ti) u(t 2 ) u(t 3 ) U(t4) v I 

R% ,R% ,R a % ,R a % ,R a % ,u(t) 

V u(ti) u(t 2 ) u(t 3 ) u(t 4 ) U(t5j v 1 



o<t<h 
h < t < t 2 
t 2 <t<t 3 
t 3 <t < u 
U<t <t 5 
t > t 5 



< 



'v(t) 

v(ti) v(t 2 ) v(t 3 ) v(t 4 ) v 1 
. R v(t\) R v{tl) R v(tl) R v(tt) K(tt) V (*) 



< t < t! 

t 1 <t<t 2 

t 2 <t<t 3 

t 3 <t< t A 
U < t < h 
t > t 5 




< t < ii 
h < t < t 2 
t 2 <t<t 3 
t 3 <t < t 4 
U < t < t 5 
t > t 5 
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fb(t) 



^S)b(t) 

R w(t\) R w{tl) R w{t 3 ) 

p — Ql D— «2 p— «3 p— Q4 
^(tO^V^) ^(43)^(44) 

p— «1 p— «2 p— «3 P _Q 4 P~"5 
^(ti ) -"V(t2 ) ^V(t 3 ) rt v(t5) 



b(t) 



b(f) 



< * < ii 
t 1 <t<t 2 
t 2 <t<t 3 
t 3 < t < U 
u < t < t 5 
t > t 5 



Some explanation is in order. Consider t > t 5 . For such a time, we first create a coun- 
terclockwise bend in the image of ^y(t) on Mi of angle about the axis passing through 
u(£ 5 ), and, simultaneously, a counterclockwise bend in the image on M 2 of angle — a^ about 
the axis passing through v(t 5 ). At the same time, we also apply these rotations to the unit 
tangent vectors a(t) and b(£), respectively These rotations create discontinuities at t 5 for 
the functions a a (£) and h a (t), but the path through M, the space of equivalence classes, 
remains continuous. We then repeat the same bending process at then at £3, then at t 2 , 
and finally at t\, creating a new absolutely continuous path 7 Q (t) through M whose images 
on Mi and M 2 are bent at five points. 

If all «j is sufficiently near 0, then 7 a (T) will be a state in M near n to use the coordinate 
system we defined above. Let Uj(a) = Wj (-y a (T)), the jth coordinate of the image of 7° at 
time T under the coordinate system defined above. The map 

uj(ai,a 2 ,a 3 ,a 4 ,a 5 ) = (u^a) , u 2 (a) , u 3 (a) , u 4 (a) , u 5 (a)) 

is clearly smooth, and if its derivative 





dui 3 ' 


dai 


dai 


duii 


dui 5 


.das 





is non-singular, then u maps open sets in IR 5 to open sets in M. As a consequence, there 
exists T > T and a G M 5 such that co(a) are the coordinates of j(T) under the coordinate 
system described above. The arc length of 7" from time to time t, however, is still t, since 
the construction of 7° above does not change arc length. Therefore 7 cannot be the length 
minimizing curve from p to ^f(T), since 7° is shorter, but this leads to a contradiction. Thus 
the derivative of u must be singular at a = (cti, a 2 , a 3 , a 4 , a 5 ) = (0, 0, 0, 0, 0) if 7 is a length 
minimizing geodesic. 

To take advantage of this observation, we examine the derivatives To compute 



dujj_ 
da; ' 



dai 

we fix all angles at except ctj. Thus computing these derivatives is reduced to 



understanding the actions of rotating counterclockwise about the axes passing through v^U) 
and v(tj). These actions will only depend on the locations of u(tj) and v(t;). Specifically, 
we compute 

dxi a {t) n ^ dv Q (t) 



da,- 



(U 2 (U), -Mi(ti),0) 



da,- 



(-rv 2 (ti),rvi(U),0) 
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da a (t) 



(0,u 3 (ti),u 2 {ti)) 



db a {t) 



(O.-Vafe),-^*)) 



Immediately this gives us derivatives for the first four coordinates: 



doti 

duj 3 
da>i 



U 2 (U 



-rv 2 (ti 



du 2 _ 
don 



-ui(U) 



rvx(ti 



The derivative of the last coordinate o>5 is only slightly more difficult: it requires just an 
application of the product rule and the fact the j a (t) passes through n at time T: 

du 5 



da. 



u 3 (U) - v 3 (U) 



Using the calculations above, we can fill in the matrix 





du>5_~ 




da i 


duji 




-da 5 


da 5 _ 



U 2 {ti) 


-Uiiti) 


-rv 2 {t{) 


rvi(ti) 




- v 3 (h 


u 2 (t 2 ) 


-Ui(t 2 ) 


-rv 2 (t 2 ) 


rv!(t 2 ) 


u 3 {t 2 ) 


- v 3 (t 2 


u 2 {t 3 ) 


~Ui{t 3 ) 


-rv 2 (t 3 ) 


rvi(t 3 ) 


u 3 {t 3 ) 


- v 3 (t 3 


u 2 {U) 


-Ml(t 4 ) 


-rv 2 (U) 


rvi(f 4 ) 


u 3 {h) 


- v 3 (U 


u 2 {h) 


-ui(h) 


-rv 2 (t 5 ) 


rvi(t 5 ) 


u 3 {t b ) 


- v 3 (t 5 



As we argued before, the above matrix must be singular. Therefore, there must exist con- 
stants ki, k 2 , k 3 , A; 4 , and k 5 , not all equal to 0, such that, for all i = 1,2,3,4,5, we have 
that 

hu^ti) + k 2 u 2 (U) + k 3 vi{ti) + hv 2 (ti) + k 5 (u 3 (ti) - v 3 (ti)) = 

The constants ki are determined by four of the times tj, and replacing the remaining time 
with any other t G [0,T] yields the same equation. Thus, for any time t G [0, T], we must 
have 

hui{t) + k 2 u 2 {t) + k 3 vi(t) + hv 2 (t) + k 5 (u 3 (t) - v 3 (t)) = 

Because of the homogeneity of M, there is nothing unique about n, so a similar plane 
equation holds for length minimizing geodesies passing through any other point in M. Thus, 



if 7 is a length minimizing geodesic, there exists a 5-plane P in 



x 



such that all points 



of contact (u(t),v(£)) G 



for 7 lie on P. 
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